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Abstract. Photosynthesis is arguably the fundamental process of Life, since it 
enables energy from the Sun to enter the food-chain on Earth. It is a remarkable non- 
equilibrium process in which photons are converted to many-body excitations which 
traverse a complex biomolecular membrane, getting captured and fueling chemical 
reactions within a reaction-center in order to produce nutrients. The precise nature of 
these dynamical processes - which lie at the interface between quantum and classical 
behaviour, and involve both noise and coordination - are still being explored. Here we 
focus on a striking recent empirical finding concerning an illumination-driven transition 
in the biomolecular membrane architecture of Rsp. Photometricum purple bacteria. 
Using stochastic realisations to describe a hopping rate model for excitation transfer, we 
show numerically and analytically that this surprising shift in preferred architectures 
can be traced to the interplay between the excitation kinetics and the reaction center 
dynamics. The net effect is that the bacteria profit from efficient metabolism at low 
illumination intensities while using dissipation to avoid an oversupply of energy at high 
illumination intensities. 
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1. Introduction 

In addition to its intrinsic interest as one of Nature's oldest and most important 
processes, photo-energy conversion is of great practical interest given Society's pressing 
need to reduce reliance on fossil fuels by exploiting alternative energy production. 
Photosynthesis maintains the planet's oxygen and carbon cycles in equilibrium [U EJ |3] 
and efficiently converts sunlight [U EJ [6], while the possibility of its in vivo study 
provides a fascinating window into the aggregate effect of millions of years of natural 
selection. Among the most widespread photosynthetic systems are purple bacteria 
Rsp. Photometricum which manage to sustain their metabolism even under dim light 
conditions within ponds, lagoons and streams [7]. They absorb light through antenna 
structures in the biomolecular Light Harvesting complex 2 (LH2), and transfer the 
electronic excitation along the membrane to Light Harvesting complexes 1 (LH1) which 
each contain a Reaction Center (RC) complex. If charge carriers are available (i.e. the 
RC is in an open state), then the resulting reactions will feed the bacterial metabolism. 

It was recently observed [1] that the photosynthetic membranes in Rsp. 
Photometricum adapt to the light intensity conditions under which they grow. 
Illuminated under High Light Intensity (HLI) (Iq w 100W/m 2 where Iq is the 
growing light intensity), membranes grow with a ratio of antenna-core complexes (i.e. 
stoichiometry) LH2/LH1^ 3.5-4. For Low Light Intensity (LLI) (J Q 10W/m 2 ), this 
ratio increases to 7-9. The features that reveal an unexpected change in the ratio of 
harvesting complexes, in bacterias grown under HLI and LLI are shown in FigUJa) 
and FigJJJb), respectively. Here we present a quantitative theory to explain this 
adaptation in terms of a dynamical interplay between excitation kinetics and reaction- 
center dynamics. In particular, the paper lays out the model, its motivation and 
implications, in a progressive manner in order to facilitate understanding. Although 
our model treats the excitation transport as a noisy, classical process, we stress that the 
underlying quantities being transported are quantum mechanical many-body excitations 
0, [9j . The membrane architecture effectively acts as a background network which loosely 
coordinates the entire process. 

2. Structure of complexes and excitation kinetics in small networks 

Figure [2] summarizes the relevant biomolecular complexes in purple bacteria Rsp. 
Photometricum [10] . together with timescales governing the excitation kinetics and 
reaction center dynamics. Each LH2 can absorb light with wavelengths 839 to 846 nm, 
while LH1 absorbs maximally at 883 nm. The LH1 forms an ellipse which completely 
surrounds the reaction center (RC) complex. Within the RC, a dimer of bacterio- 
chlorophyls (BChls) known as the special pair P, can be excited. The excitation (P*) 
induces ionization (P + ) of the special pair, and hence metabolism. The initial photon 
absorption is proportional to the complex cross-sections, which have been calculated 
for LH1 and LH2 complexes [IT] . With n(X) incident photons of wavelength A, an 
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Figure 1. Empirical architectures for (a) HLI and (b) LLI membranes, displaying 
LH2s (small blue circles), LHls (big green circles) and RCs (orange dots). 




Figure 2. Schematic of the biomolecular photosynthetic machinery in purple bacteria, 
together with relevant inter-complex mean transfer times ty, dissipation rate jd, and 
normalized light intensity rate 71(2) 

18 W/m 2 light intensity yields a photon absorption rate for circular LH1 complexes 
in Rb. sphaeroides [12] given by J n(A)<7LHi(A)<iA = 18s -1 , where <7 L hi is the LH1 
absorption cross section. For LH2 complexes, the corresponding photon capture rate is 
10s -1 . Extension to other intensity regimes is straightforward, by normalizing to unit 
light intensity. The rate of photon absorption normalized to 1 W/m 2 intensity, will 
be 71 = y§ = Is" 1 for an individual LH1, and 72 = y§ = 0.55s -1 for individual LH2 
complexes. The complete vesicle containing several hundreds of complexes, will have an 
absorption rate 7^ = 1(71 A^i + 72./V2) where N±(2) is the number of LH1 (LH2) complexes 
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in the vesicle, and I is the light intensity. The number of RC complexes is therefore 
also equal to N\. Excitation transfer occurs through induced dipole transfer, among 
BChls singlet transitions. The common inter-complex BChl distances 20-100 A [TJ [13] 
cause excitation transfer to arise through the Coulomb interaction on the picosecond 
time-scale [H], while vibrational dephasing destroys coherences within a few hundred 
femtoseconds [l5j[T6]. The Coulomb interaction de-excites an initially excited electron in 
the donor complex while simultaneously exciting an electron in the acceptor complex. 
As dephasing occurs, the donor and acceptor phase become uncorrelated. Transfer 
rate measures from pump-probe experiments agree with generalized Forster calculated 
rates [14] . assuming intra-complex derealization. LH2— >LH2 transfer has not been 
measured experimentally, although an estimate of t 22 = 10 ps has been calculated [14] . 
LH2— > LH1 transfer has been measured for R. Sphaeroides as t 2 i = 3.3ps [17]. Due to 
formation of excitonic states [HI E], back-transfer LH1— > LH2 is enhanced as compared 
to the canonical equilibrium rate for a two-level system, up to a value of t i2 = 15.5ps. 
The LH1— 7-LH1 mean transfer time tu has not been measured, but generalized Forster 
calculation [18] has reported an estimated mean time t u of 20 ps. LH1— > RC transfer 
occurs due to ring symmetry breaking through second and third lowest exciton lying 
states [19], as suggested by agreement with the experimental transfer time of 35-37 ps 
at 77 K [201 HI] • Increased spectral overlap at room temperature improves the transfer 
time to ti t Rc = 25 ps as proposed by [22J. A photo-protective design makes the back- 
transfer from an RC's fully populated lowest exciton state to higher-lying LH1 states 
occur in a calculated time of £_rc*,i =8.1 ps [TH], close to the experimentally measured 
7-9 ps estimated from decay kinetics after RC excitation |23j . The first electron transfer 
step P* — > P + occurs in the RC within t+ =3 ps, used for quinol (QbH 2 ) production 
[14] . Fluorescence, inter-system crossing, internal conversion and further dissipation 
mechanisms, have been included within an effective single lifetime I/70 of 1 ns [18J. 
Due to the small absorption rates in 7^, two excitations will only rarely occupy a single 
harvesting structure - hence it is sufficient to include the ground s = and one exciton 
states s = 1 for each harvesting complex. 

We now introduce the theoretical framework that we use to describe the excitation 
transfer, built around the experimental and theoretical parameters just outlined. In 
the first part of the paper, our calculations are all numerical - however we turn to an 
analytic treatment in the latter part of the paper. We start by considering a collective 
state with N = N 2 + 2Ni sites - resulting from JVi LHls, N 2 LH2s and hence N x RC 
complexes in the vesicle - in terms of a set of states having the form {si,..,Sat} i n 
which any complex can be excited or unexcited, and a maximum of N excitations can 
exist in the membrane. If only excitation kinetics are of interest, and only two states 
(i.e. excited and unexcited) per complex are assumed, the set of possible states has 2^ 
elements. We introduce a vector p = (pi, p 2 , P2 N ) in which each element describes the 
probability of occupation of a collective state comprising several excitations. Its time 
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M) = ^#^777 ■ (2) 



evolution obeys a master equation 

niV 

d tPl (t) = Y.G h3 p 3 (t). (1) 

Here Gij is the transition rate from a site % to a site j. Since the transfer rates do 
not depend on time, this yields a formal solution p{t) = e Gt p(0). Small absorption 
rates lead to single excitation dynamics in the whole membrane, reducing the size of 
p(t) to the total number of sites N. The probability to have one excitation at a given 
complex initially, is proportional to its absorption cross section, and can be written as 
p(0) = -^(71, •••,72, 0, ..), where subsets correspond to the N\ LHls, the N 2 LH2s and 

JVi N 2 Ni 

the Ni RCs respectively. Our interest lies in p^ which is the normalized probability to 
find an excitation at a complex, given that at least one excitation resides in the network: 

Pk{t) 

In order to appreciate the effects that network architecture might have on the 
model's dynamics, we start our analysis by studying different arrangements of complexes 
in small model networks, focusing on architectures which have the same amount of 
LH1, LH2 and RCs as shown in the top panel of Figj3](a), (b) and (c). The bottom 
panel Figj3] (d)-(f ) shows that p^ values for RC, LH1 and LH2 complexes, respectively. 
FigJHJ^d) shows that the highest RC population is obtained in configuration (c), followed 
by configuration (a) and (b) whose ordering relies in the connectedness of LHls to 
antenna complexes. Clustering of LHls will limit the number of links to LH2 complexes, 
and reduce the probability of RC ionization. For completeness, the probability of 
occupation in LH1 and LH2 complexes (Figs|3fe) and (f), respectively), shows that 
increased RC occupation benefits from population imbalance between LH1 enhancement 
and LH2 reduction. As connections among antenna complexes become more favored, 
the probability of finding an excitation on antenna complexes will become smaller, while 
the probability of finding excitations in RCs is enhanced. 

This discussion of simple network architectures, provides us with a simple 
platform for testing the notion of energy funneling, which is a phenomenon that is 
commonly claimed to arise in such photosynthetic structures. We start with a minimal 
configuration corresponding to a basic photosynthetic unit: one LH2, one LH1 and 
its RC. Figure H^a) shows that excitations will mostly be found in the LH1 complex, 
followed by occurrences at the LH2 and lastly at the RC. Figure 0Jb) shows clearly 
the different excitation kinetics which arise when the RC is initially unable to start the 
electron transfer P* — > P + , and then after 15ps the RC population increases with 
respect to the LH2's. This confirms that the energy funneling concept is valid for these 
small networks [El HE], i.e. excitations have a preference to visit the RC (ti jRC = 25ps) 
as compared to being transferred to the light-harvesting complexes (tu = 15.5ps). 
However, in natural scenarios involving entire chromatophores with many complexes, we 
will show that energy funneling is not as important due to increased number of available 
states, provided from all LH2s surrounding a core complex. 
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Figure 3. Top panel: Three example small network architectures. The bottom panel 
shows the normalized probabilities for finding an excitation at an RC (see (d)), an LH1 
(see (e)), or an LH2 (see (f)). In panels (d)-(f), we represent these architectures as 
follows: (a) is a continuous line; (b) is a dotted line; (c) is a dashed line. 
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Figure 4. Normalized probabilities pfc for finding the excitation at an LH2 (dashed), 
LH1 (dotted) or at an RC (continuous), for (a) t+ = 3ps, and (b) t + — > oo. Crosses 
are the results from the Monte Carlo simulation. 



Given the large state-space associated with such multiple complexes, our subsequent 
model analysis will be based on a discrete-time random walk for excitation hopping 
between neighboring complexes. In particular, we use a Monte Carlo method to 
simulate the events of excitation transfer, the photon absorption and dissipation, and 
the RC electron transfer. We have checked that our Monte Carlo simulations accurately 
reproduce the results of the population-based calculations described above, as can be 
seen from FigsJUa) and (b). 
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Figure 5. Contour plots for dissipation in LLI (a) and HLI (b) membranes. Greater- 
contrast means higher values for dissipation. The simulation is shown after 10 
excitations were absorbed by the membrane with rate ja- 

3. Performance measures of complete chromatophore vesicles 

We now turn to discuss the application of the model to the empirical biological structures 
of interest, built from the three types of complex k (LH1, k=l; LH2, k=2; RC, k=3). 
In particular, we have carried out extensive simulations to investigate the role of the 
following quantities in the complete chromatophore vesicles: 

• Adjacency geometry of LHls and LH2s. The LH2s are more abundant than LHls 
and both complexes tend to form clusters, while LH2s are also generally found 
surrounding the LHls. 

• The average time an excitation spends in complex type k. 

• The probability pn k of finding an excitation in complex type k. 

• Dissipation di which measures the probability for excitations to dissipate at site i, 
from which the probability of dissipation in core or antenna complexes can be 
obtained by adding all di concerning complex type k. 

• The sum over all complexes of the dissipation probability, which gives the 
probability for an excitation to be dissipated. The efficiency of the membrane 
is the probability of using an excitation in any RC, i.e. 77 = 1 — J2i di- 

Figure E^a) shows that the membrane grown under low light intensity (LLI) has 
highly dissipative clusters of LH2s, in contrast to the uniform dissipation in the high 
light intensity (HLI) membrane (see Figj5](b)). This result is supported by a tendency 
for excitations to reside longer in LH2 complexes far from core centers (not shown), 
justifying the view of LH2 clusters as excitation reservoirs. However, for LLI and 
HLI, the dissipation in LH1 complexes is undistinguishable. In Table [1] we show the 
observables obtained using our numerical simulations. These show that: 
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(i) Funneling of excitations: 

• The widely held view of the funneling of excitations to LH1 complexes, turns 
out to be a small network effect, which by no means reflects the behavior 
over the complete chromatophore. Instead, we find that excitations are found 
residing mostly in LH2 complexes. 

• Since a few LH2s surround each LH1, the mean residence times t/~ in all 
complexes is very similar. 

(ii) Dissipation and performance: 

• Excitations are dissipated more efficiently in individual LH1 complexes, since 

Hi v. £2 
Ni N2 ' 

• Dissipation in a given complex type depends primarily on its relative 
abundance, since jjr ~ jfr 

• HLI membranes are more efficient than LLI membranes. 
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Table 1. Residence time tk (in picoseconds), dissipation D^, residence probability 
Pn k: unitary dissipation per complex (xl0~ 3 ), on k = {1,2} corresponding to N\ 
LH1 and N2 LH2 complexes respectively. Stoichiometry s and efficiency rj are also 
shown. 



For the present discussion, the most important finding from our simulations is 
that the adaptation of purple bacteria does not lie in the single excitation kinetics. 
In particular, LLI membranes are seen to reduce their efficiency globally at the point 
where photons are becoming scarcer - hence the answer to adaptation must lie in some 
more fundamental trade-off (as we will later show explicitly). Due to the dissimilar 
timescales between millisecond absorption [12] and nanosecond dissipation [H] , multiple 
excitation dynamics are also unlikely to occur within a membrane. However we note 
that simulations involving multiple excitations, that include blockade (FigH^a)) in which 
two excitations can not occupy the same site, does not appreciably lower the efficiency 77 
up to thirty excitations. We find that annihilation (FigJH^b)), in which two excitations 
annihilate when they occupy the site at the same time, diminishes the membrane's 
performance equally in both HLI and LLI membranes. 

Our findings above show that the explanation for the observed architecture 
adaptations (HLI and LLI) neither lies in the frequently quoted side-effect of multiple 
excitations, nor in the excitation dynamics alone. Instead, as we now explain, the answer 
as to how adaptation can prefer the empirically observed HLI and LLI structures under 
different illumination conditions, lies in the interplay between the excitation kinetics and 
reaction-center cycling dynamics. By virtue of quinones-quinol and cytochrome charge 
carriers, the RC dynamics features a 'dead' (or equivalently 'busy') time interval during 
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Figure 6. Efficiency of multiple excitation dynamics: (a) blockade and (b) 
annihilation mechanisms, for LLI (crosses) and HLI (diamonds) membranes, no 
corresponds to the initial number of excitations in each realization. 




Figure 7. Reaction Center cycle, showing double reduction of the special pair P 
together with formation of quinol QbH 2 - There is a dead time r on the millisecond 
time-scale, before a new quinone Qb becomes available. 



which quinol is produced, removed and then a new quinone becomes available [2U |25] . 
A single oxidation P* — > P + will produce QbH in the reaction Qb — >■ Qb ~^ QbH, and 
a second oxidation will produce quinol QbHi in the reaction QbH — > QbH~ — > QbH 2 . 
Once quinol is produced, it leaves the RC and a new quinone becomes attached. The 
cycle is depicted in Fig. [7J and is described in the simulation algorithm by closing an RC 
for a time r after two excitations form quinol. This RC cycling time r implies that at any 
given time, not all RCs are available for turning the electronic excitation into a useful 
charge separation. Therefore, the number of useful RCs decreases with increasing r. Too 
many excitations will rapidly close RCs, implying that any subsequently available nearby 
excitation will tend to wander along the membrane and eventually be dissipated - hence 
reducing rj. For the configurations resembling the empirical architectures (FigJTJ, this 
effect is shown as a function of r in Fig. [8(a) yielding a wide range of RC-cycling times 
at which LLI membrane is more efficient than HLI. Interestingly, this range corresponds 
to the measured time-scale for r of milliseconds [2U [25] , and supports the suggestion 
that bacteria improve their performance in LLI conditions by enhancing quinone-quinol 
charge carrier dynamics as opposed to manipulating exciton transfer. A recent proposal 
[27] has shown numerically that the formation of LH2 para-crystalline domains produces 
a clustering trend of LH1 complexes with enhanced quinone availability - a fact that 
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would reduce the RC cycling time. However, the crossover of efficiency at r m 3 ms 
implies that even if no enhanced RC-cycling occurs, the HLI will be less efficient than the 
LLI membranes on the observed r time-scale. The explanation is quantitatively related 
to the number N Q of open RCs. Figs. WJo), (c) and (d) present the distribution p(N ) 
of open RCs, for both HLI and LLI membranes and for the times shown with arrows in 
FigJH^a). When the RC-cycling is of no importance (Fig. [S(b)) almost all RCs remain 
open, thereby making the HLI membrane more efficient than LLI since having more 
(open) RCs induces a higher probability for special pair oxidation. Near the crossover 
in Fig. [8], both membranes have distributions p{N Q ) centered around the same value 
(Fig. El(c)), indicating that although more RCs are present in HLI vesicles, they are 
more frequently closed due to the ten fold light intensity difference, as compared to LLI 
conditions. Higher values of r (Fig. 0(d)) present distributions where the LLI has more 
open RCs, in order to yield a better performance when photons are scarcer. Note that 
distributions become wider when RC cycling is increased, reflecting the mean-variance 
correspondance of Poissonian statistics used for simulation of r. Therefore the trade- 
off between RC-cycling, the actual number of RCs and the light intensity, determines 
the number of open RCs and hence the performance of a given photosynthetic vesicle 
architecture (i.e. HLI versus LLI). Guided by the Monte Carlo numerical results, we 
develop in Sec. [5] an analytical model (continuous lines in FigJH]) that supports this 
discussion. 

For completeness, we now quantify the effect of incident light intensity variations 
relative to the light intensity during growth, with both membranes having r = 3ms. 
The externally applied light intensity I/Io, which corresponds to the ratio between the 
actual (J) and growth (Jo) light intensities, is varied in Fig. [9](a). The LLI membrane 
performance starts to diminish well beyond the growth light intensity, while the HLI 
adaptation starts diminishing just above Iq due to increased dissipation. The crossover 
in efficiency at / ~ Iq results from the quite different behaviors of the membranes 
as the light intensity increases. In particular, in LLI membranes excess photons are 
readily used for bacterial metabolism, and HLI membranes exploit dissipation in order 
to limit the number of processed excitations. Figs. Mjo), (c) and (d) verify that 
performance of membranes heavily depends on the number of open RCs. For instance, 
membranes subject to low excitation intensity (Fig. E](b)) behave similarly to that 
expected for fast RC cycling times (Fig. [S^a)). The complete distributions, both for 
HLI and LLI conditions, shift to lower N a with increased intensity in the same manner 
as that observed with r. Even though these adaptations show such distinct features in 
the experimentally relevant regimes for the RC-cycling time and illumination intensity 
magnitude [H[2ll[25], FigsJB^c) and (d) show that the distributions of open RCs actually 
overlap. Despite the fact that the adaptations arise under different environmental 
conditions, the resulting dynamics of the membranes are quite similar. Note that within 
this parameter subspace of I and r, the LLI membrane may have a larger number of 
open RCs than the HLI adaptation. In such a case, the LLI membrane will perform 
better than HLI with respect to RC ionization. The inclusion of RC dynamics implies 
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Figure 8. (a) Monte Carlo calculation of efficiency 77 of HL1 (diamonds) and LLI 
(crosses) grown membranes, as a function of the RC-cycling time r. Continuous lines 
give the result of the analytical model, (b), (c) and (d) show the distributions p(N ) 
of the number of open RCs for the times shown with arrows in the main plot for HLI 
(filled bars) and LLI (white bars). 



that the absorbed excitation will not find all RCs available. Instead, a given amount of 
closed RCs will eventually alter the excitation's fate since probable states of oxidization 
are readily reduced. In a given lifetime, an excitation will find (depending on r and 
I) a number of available RCs - which we refer to as effective stoichiometry - which is 
different from the actual number reported by Atomic Force Microscopy [TJ [13] . 



4. Clustering trends 

The empirical AFM investigations show two main features which highlight the 
architectural change in the membrane as a result of purple bacteria's adaptation: the 
stoichiometry variation and the trend in clustering. Fig. [10] shows the importance of 
the arrangement of complexes by comparing architectures (b), (c), (d) and (e), all of 
which have a stoichiometry which is consistent with LLI vesicles. Fig. [IDT a) shows the 
difference between a given membrane's efficiency 77, and the mean of all the membranes 
77, i.e. A77 = 77 — 77. The more clustered the RCs, the lower the efficiency in the short 
r domain. As RC cycling is increased, (b) becomes the least efficient while all other 
configurations perform almost equally. The explanation comes from the importance 
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Figure 9. (a) Results from our Monte Carlo calculation of efficiency r\ of HLI 
(diamonds, 7 = 100W/m 2 ) and LLI (crosses, Jo = 10W/m 2 ) membranes, as a function 
of incident light intensity J/Jo- Continuous lines give the result of the analytical model. 
Panels (b), (c) and (d) show the distribution p{N ) of the number of open RCs for 
light intensities corresponding to arrows in the main plot. HLI are shaded bars and 
LLI are white bars. 

of the number of open RCs: As r gets larger, many RCs will close and the situation 
becomes critical at r « 3ms, where r] decreases rapidly. Configurations (c), (d) and (e) 
all have the same number of RCs (i.e. 44), and the distribution of open RCs is almost 
the same in each case for any fixed RC cycling time. By contrast, (b) has fewer RCs 
(i.e. 36). Therefore when r is small, sparser RCs and exciton kinetics imply that the 
membrane architecture (b) will have better efficiency than (c) and (e). The effect of 
the arrangement itself is lost due to slower RC dynamics, and the figure of merit that 
determines efficiency is the number of open RCs, which is lower for (b). 

To summarize so far, we find that the arrangement of complexes changes slightly 
the efficiency of the membranes when no RC dynamics is included - but with RC 
dynamics, the most important feature is the number of open RCs which is smaller for 
(b). The nearly equal efficiency over the millisecond r domain, emphasizes the relative 
insensitivity to the complexes' geometrical arrangement. The slower the RC cycling, the 
more evenly available RCs will be dispersed in clustered configurations, resembling the 
behavior of sparse RC membranes. Incoming excitations in clustered configurations may 
quickly reach a cluster bordering closed RCs, but must then explore further in order to 
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Figure 10. (a) A-q is presented for the following membrane configurations: (b) Rsp. 
Photometricum bacteria; (c) Rs. Palustris bacteria; (d) completely unclustered vesicle; 
(e) fully clustered vesicle. The symbols are crosses, circles, diamonds and boxes, 
respectively 



generate a charge separation. Although the longer RC closing times make membranes 
more prone to dissipation and decreased efficiency, it also makes the architecture less 
relevant for the overall dynamics. The relevant network architecture instead becomes the 
dynamical one including sparse open RCs, not the static geometrical one involving the 
actual clustered RCs. The inner RCs in clusters are able to accept excitations as cycling 
times increase, and hence the RCs overall are used more evenly. This implies that there 
is little effect of the actual configuration, and explains the closeness of efficiencies for 
different arrangements in the millisecond range. 

5. Global membrane model incorporating excitation kinetics and 
RC-cycling 

Within a typical fluorescence lifetime of 1 ns, a single excitation has travelled hundreds 
of sites and explored the available RCs globally. The actual arrangement or architecture 
of the complexes seems not to influence the excitation's fate, since the light intensity 
and RC cycling determine the number of open RCs and the availability for P oxidation. 
This implies that the full numerical analysis of the excitation kinetics, while technically 
more accurate, may represent excessive computational effort - either due to the size of 
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the state space within the master equation approach, or the number of runs required for 
ensemble averages with the stochastic method. In addition, within neither numerical 
approach is it possible to deduce the direct functional dependence of the efficiency 
on the parameters describing the underlying processes. To address these issues, we 
present here an alternative rate model which is inspired by the findings of the numerical 
simulations, but which (1) globally describes the excitation dynamics and RC cycling, 
(2) leads to analytical expressions for the efficiency of the membrane and the rate of 
quinol production, and (3) sheds light on the trade-off between RC-cycling and exciton 
dynamics [26] . 

We start with the observation that absorbed excitations are transferred to RCs, 
and finally ionize special pairs or are dissipated. At any given time Ne excitations will 
be present in the membrane. The rate at which they are absorbed is 7,4. Excitations 
reduce quinone Qb in the membrane at RCs due to P oxidation at a rate XcNe, or 
dissipate at a rate JdNe- Both processes imply that excitations leave the membrane at 
a rate ^f- = —XcNe —^dNe- Here Ac is the inverse of the mean time trc at which an 
excitation yields a charge separation at RCs when starting from any given complex, and 
it depends on the current number of open RCs given by iV . The RC cycling dynamics 
depend on the rate at which RCs close, i.e. —^^-Ne where the 1/2 factor accounts 
for the need for two excitations to produce quinol and close the RC. The RCs open at 
a rate 1/r, proportional to the current number of closed RCs given by iV\ — N a . Hence 
the RC-excitation dynamics can be represented by two nonlinear coupled differential 
equations: 

^f=- (A c (iV ) + lD )N E + 1a (3) 

In the stationary state, the number of absorbed excitations in a time interval At and 
the number of excitations used to produce quinol urc, are given by 

n A = 7 A At (5) 

n RC = X c (N )N E At (6) 

yielding an expression for the steady-state efficiency r\ = n RC l n A- 

Xc(N )N E _ 
r, = . (7) 

7a 

These equations can be solved in the stationary state for Ne and N Q , algebraically 
or numerically, if the functional dependence of Xc{N a ) is given. It is zero when 
all RCs are closed, and a maximum X c when all are open. Making the functional 
dependence explicit, Xc{N ), Fig. UTT a) presents the relevant functional form for HLI 
and LLI membranes together with a linear and a quadratic fit. The dependence on the 
rate of quinone reduction Xc(N ) requires quantification of the number of open RCs, 
with a notation where the fitting parameter comprises duplets where first and second 
components relate to the HLI and LLI membranes being studied. Figure [TTT a) shows 
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Figure 11. (a) Numerical results showing the rate of ionization Xc(N ) of an RC 
for HLI (diamonds) and LLI (crosses) membranes, together with a quadratic (dashed 
line) and linear (continuous) dependence on the number of closed RCs (Ni —N Q ). The 
fitting parameters for a + bN Q are a = {15.16, 7.72}ns _1 , b = {—0.21, —0.21} ns _1 ; and 
for a + bN + cNl, a = {16.61, 8.211ns- 1 , b = {-0.35, -0.33}, and c = {3.6, 1.5}/is"\ 
for HLI and LLI membranes respectively, (b) rj as function of r and a = I /Iq, obtained 
from the complete analytical solution for LLI (white) and HLI (grey) membranes 

that \c{N ) favors a quadratic dependence of the form Ac(iV ) = a + bN a + cN%. A 
linear fit 

A c (iV ) = A° c (^) (8) 

smears out the apparent power-law behavior with fit value Ni = {70.72, 35.71}, in close 
agreement with the HLI and LLI membranes which have 67 and 36 RCs respectively. 
The linear fit Ac(iV ) can used to generate an analytical expression for rj as follows: 

V(r, 7 aC0) = -4t" { 2N ^c + Id) + lA\° c r- (9) 

yj4N?(X° c + Id) 2 + ±N llA \ Q c { lD - X° c )r + (7aA° c t) 2 | . 

We found no analytical solution for rj in the case where Xc(N ) has a power law 
dependence. In the limit of fast RC cycling-time (r— >0), r\ has the simple form 
j) — (1 + 7£)/A^) _1 . If all transfer paths are summarized by A^, this solution illustrates 
that rj > 0.9 if the transfer-P reduction time is less than a tenth of the dissipation 
time, not including RC cycling. As can be seen in Figs. |8] and the analytical 
solution is in good quantitative agreement with the numerical stochastic simulation, 
and provides support for the assumptions that we have made. Moreover, our theory 
shows directly that the efficiency is driven by the interplay between the RC cycling time 
and light intensity. Figure fTlT b) shows up an entire region of parameter space where LLI 
membranes are better than HLI in terms of causing P ionization, even though the actual 
number of RCs that they have is smaller. In view of these results, it is interesting to 
note how clever Nature has been in tinkering with the efficiency of LLI vesicles and the 
dissipative behavior of HLI adaptation, in order to meet the needs of bacteria subject 
to the illumination conditions of the growing environment. 
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6. Bacterial metabolic demands: A stoichiometry prediction from our 
model 

Photosynthetic membranes must provide enough energy to fulfill the metabolic 
requirements of the living bacteria. In order to quantify the quinol output of the vesicle, 
we calculate the quinol rate 

which depends directly on the excitations that ionize RCs tirc- The factor | accounts for 
the requirement of two ionizations to form a single quinol molecule. Fig. IT27 a) shows the 
quinol rate as a function of RC cycling time, when membranes are illuminated in their 
respective grown conditions. If RC cycling is not included (r — > 0) the tenfold quinol 
output difference suggests that the HLI membrane could increase the cytoplasmic pH to 
dangerously high levels, or that the LLI membrane could starve the bacteria. However, 
the bacteria manage to survive in both these conditions - and below we explain why. 

In the regime of millisecond RC cycling, the quinol rate in HLI conditions decreases 
which is explained by dissipation enhancement acquired from only very few open RCs. 
Such behavior in LLI conditions appears only after several tens of milliseconds. The fact 
that no crossover occurs in quinol rate for these two membranes, suggests that different 
cycling times generate this effect. The arrows in Fig. [I~27 a) correspond to times where 
a similar quinol rate is produced in both membranes, in complete accordance with 
numerical studies where enhanced quinone diffusion lessens RC cycling times [27] in 
LLI adaptation. Although these membranes were grown under continuous illumination, 
the adaptations themselves are a product of millions of years of evolution. Using RC 
cycling times that preserve quinol rate in both adaptations, different behaviors emerge 
when the illumination intensity is varied (see Fig. [T2|) . The increased illumination is 
readily used by the LLI adaptation, in order to profit from excess excitations in an 
otherwise low productivity regime. On the other hand, the HLI membrane maintains 
the quinone rate constant, thereby avoiding the risk of pH imbalance in the event that 
the light intensity suddenly increased. We stress that the number of RCs synthesized 
does not directly reflect the number of available states of ionization in the membrane. 
LLI synthesizes a small amount of RCs in order to enhance quinone diffusion, such that 
excess light intensity is utilized by the majority of special pairs. In HLI, the synthesis 
of more LH1-RC complexes slows down RC-cycling, which ensures that many of these 
RCs are unavailable and hence a steady quinol supply is maintained independent of 
any excitation increase. The very good agreement between our analytic results and the 
stochastic simulations, yields additional physical insight concerning the stoichiometries 
found experimentally in Rsp. Photometricum\\\. In particular, the vesicles studied 
repeatedly exhibit the same stoichiometries, s ~ 4 for HLI, and s ~ 8 for LLI 
membranes. Interestingly, neither smaller nor intermediate values are found. 

We now derive an approximate expression for the quinol production rate W 
in terms of the environmental growth conditions and the responsiveness of purple 
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Figure 12. (a) Quinol rate W in HLI (diamonds, Iq = 100W/m 2 ) and LLI (crosses, 
Iq = 10W/m 2 ) grown membranes, as a function of RC cycling time. The times shown 
with arrows are used in (b), where W is presented as a function of incident intensity. 



bacteria through stoichiometry adaptation. Following Refs. [U [27J [28], the area A Q 
of chromatophores in different light intensities can be assumed comparable. Initially 
absorption occurs with Ni LH1 complexes of area A\, and N 2 LH2 complexes of area 
A 2 , which fill a fraction / of the total vesicle area / = {A\N\ + A 2 N 2 )/Aq. This 
surface occupancy can be rearranged in terms of the number of RCs Ni, and the 
stoichiometry Ni = N 2 /s, yielding the expression / = Ni(A\ + sA 2 )/A . The fraction / 
has been shown [27] to vary among adaptations f(s), since LLI have a greater occupancy 
than HLI membranes due to para-crystalline LH2 domains in LLI. Accordingly, the 
absorption rate can be cast as 7 a = /(71 + 372) a^+^a 2 • Following photon absorption, 
the quinol production rate W = XqNe/2 depends on the number of excitations within 
the membrane in the stationary state, and on the details of transfer through the rate Ac- 
The assumed linear dependence X C (N ) requires knowledge of Xq and the stationary- 
state solution for the mean number of closed RCs through N a = N\ — 2 (^o+x c ) T ' ^~^ e 
stationary state in Eq. [3] yields Ne = ^+\ c sucn that the mean number of closed RCs 
is simply N Q = Ni — Wr. The rate A^ is the rate at which excitations oxidize any 
special pair when all RCs are open. The time to reach an RC essentially depends on 
the number of RCs and hence the stoichiometry s. \° c must be zero when no RCs are 
present (s— >oo) and takes a given value (to) -1 when the membrane is made of only LHls 
(s=0). Also the RC cycling time r is expected to vary somewhat with adaptations due 
to quinone diffusion [27], which is supported in our analysis by the condition of bounded 
metabolic demands as presented in Fig. [T2T a). The linear assumed dependence of Eq. 
|8] is cast explicitly as a function of W and s with the number of open RCs 

From Ne in the steady state, we have 

2{\c(s,W)+ lD )' 1 ' 
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which can be solved for W(s,I) 

„ w , n 1a(sJ) 1 (-. 1d\ I, 7a{s,I) 1 ^ , 7dV 2 , 7A(s,I) no , 

2W(s ' I] = ^^ + W) + W) V + a° J ] + ^T (13) 

where B(s) = ^j?~rl[~^ ■ To determine r(s), we employ a linear interpolation using 
the values highlighted by arrows in Fig. IT27 b). The requirements on Xq(s) are 
fulfilled by a fit A^(s) = (s/a + (to)) -1 , which is supported by our calculations for 
configurations of different stoichiometries (see Fig. [131( a)). The filling fraction fraction 
f(s) is assumed linear according to the experimentally found values for HLI (/ rj 0.75) 
and LLI (/ w 0.85). The resulting expression is presented in Fig. [TBT b). In the 
high stoichiometry/high intensity regime, the high excitation number would dangerously 
increase the cytoplasmic pH (HI [121 EH]- The contours shown in Fig. [TBT c) of constant 
quinol production rate W, show that only in a very small intensity range will bacteria 
adapt with stoichiometries which are different from those experimentally observed in 
Rsp. Photometricum (s ~4 and s ~8). As emphasized in Ref. [I], membranes with 
s=6 or s=2 were not observed, which is consistent with our model. More generally, our 
results predict a great sensitivity of stoichiometry ratios for 30-40 W/m 2 , below which 
membranes rapidly build up the number of antenna LH2 complexes. Very recently 
[29J, membranes were grown with 30W/m 2 and an experimental stoichiometry of 4.8 
was found. The contour of 2200 s -1 predicts a value for stoichiometry of 4.72 at such 
light intensities. This agreement is quite remarkable, since a simple linear model would 
wrongly predict s = 7.1. Our theory's full range of predicted behaviors as a function of 
light-intensity and stoichiometry, awaits future experimental verification. 

7. Concluding remarks 

We have shown that excitation dynamics alone cannot explain the empirically observed 
adaptation of light-harvesting membranes. Instead, we have presented a quantitative 
model which strongly suggests that chromatic adaptation results from the interplay 
between excitation kinetics and RC charge carrier quinone-quinol dynamics. Specifically, 
the trade-off between light intensity and RC cycling dynamics induces LLI adaptation 
in order to efficiently promote P oxidation due to the high amount of open RCs. By 
contrast, the HLI membrane remains less efficient in order to provide the bacteria with 
a relatively steady metabolic nutrient supply. 

This successful demonstration of the interplay between excitation transfer and RC 
trapping, highlights the important middle ground which photosynthesis seems to occupy 
between the fast dynamical excitation regime in which quantum-mechanical calculations 
are most relevant [H [HI [30] , and the purely classical regime of the millisecond timescale 
bottleneck in complete membranes. We hope our work will encourage further study 
of the implications for photosynthesis of this fascinating transition regime between 
quantum and classical behaviors [31]. On a more practical level, we hope that our study 
may help guide the design of more efficient solar micropanels mimicking natural designs. 
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Figure 13. (a) Dependance of Ag. on stoichiometry s of the membrane, (b) W(s, I) 
as a function of stoichiometry s variation and illumination intensity, (c) Quinol rate 
contours of W = {1900, 2000, 2100, 2200} s" 1 in black, blue, red and pink, respectively 
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Appendix 

Residence time tg h In the stochastic simulations, is the residence time of an 
excitation in the a th realization at complex i: 
V t a 

tk = — (-1J 

l^iek n Vi 

where ny i is the number of times complex % has been visited in all the stochastic 
realizations. 



Dissipation di The dissipation di measures the probability for excitations to dissipate 
at site i, and can be obtained formally from 

* = (.2) 

n A 
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where n L ~, i are the number of excitations dissipated at site i and ua is the total number 
of absorbed excitations. 



Residence probability pn k In practice within the stochastic simulations, for a given 
realization a, an excitation will be a time t\ in LHls, a time t% in LH2s and a time t% 
in RCs. The residence probability will become 

V t a 

PRk = (- 3 ) 

The total time of all realizations in complex type k is the numerator, while the 
denominator stands for the total time during which the excitations were within the 
membrane. 
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